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Abstract 

One of the many surprising results found in the mechanics of rotating systems is 
the stabilization of a particle in a rapidly rotating planar saddle potential. Besides 
the counterintuitive stabilization, an unexpected precessional motion is observed. 

In this note we show that this precession is due to a Coriolis-like force caused by 
the rotation of the potential. To our knowledge this is the first example where 
such force arises in an inertial reference frame. We also propose an idea of a simple 
mechanical demonstration of this effect. 

Introduction. 

According to Earnshaw’s theorem an electrostatic potential cannot have stable equi¬ 
libria, i.e. minima, since such potentials are harmonic functions. The theorem does 
not apply, however, if the potential depends on time; in fact, the 1989 Nobel Prize 
in physics was awarded to W. Paul [T] for his invention of the trap for suspending 
charged particles in an oscillating electric field. Paul’s idea was to stabilize the saddle 
by “vibrating” the electrostatic field, by analogy with the so-called Stephenson-Kapitsa 
pendulum m El 111 016] in which the upside-down equilibrium is stabilized by vibration 
of the pivot. Instead of vibration, the saddle can also be stabilized by rotation of the 
potential (in two dimensions); this has been known for nearly a century: as early as 
1918, L.E.J. Brouwer (1881-1966), one of the authors of the fixed point theorem in 
topology, considered stability of a heavy particle on a rotating slippery surface, mn 
Brouwer derived the equations of motion in [ 710 ; the derivation took over 3 pages. He 
then linearized the equations by discarding quadratic and higher order terms in position 
and velocity. The resulting linearized equations in the stationary frame are dans] 

X + X cos 2u:t + y sin 2ujt = 0, 

( 1 ) 

y + X sin 2ojt — y cos 2ujt = 0, 

where t is dimensionless time related to the dimensional time t via t = tD/Do, Dq = 
and where oj is the dimensionless angular velocity given by w = D/Dq- 
Equations (0) look even nicer in vector form: 

X -|- S(a;f)x = 0, (2) 

^Note that Pq is the frequency of small oscillations of the particle along the x-direction on the 
non-rotating saddle. Thus our dimensionless t measnres time in the units of the period of the above 
mentioned oscillations. 
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Figure 1: A unit mass on a saddle |71|8l|9|. 


where 


S{ujt) 


cos 2ujt sin 2ujt 
sin 2ojt — cos 2ojt 


(3) 


We note that these equations are written in the stationary frame; Brouwer actually 
derived the corresponding equations in a frame rotating with the saddle ([ZIEIE]). 
The equations ([T|), when written in the rotating frame, acquire the Coriolis and the 
centrifugal forces, but lose time-dependence since the saddle appears stationary in the 
co-rotating frame; the equivalence of the autonomous linear equations of Brouwer and 
the time-periodic equations ([T]) is well-known, see e.g. naiiiiiis]. 

We now describe another context in which ([T|) appear. 


A particle in a rotating potential. The same equations ([U - ([2]) arise also in a 
different context: they govern the motion of a unit mass in the plane under the influence 
of a potential force given by the rotating saddle potential - namely, the potential whose 
graph is obtained by rotating the graph of z = — around the z-axis with angular 

velocity oj. This planar problem is related to, but is different from that of Figure [TJ 
Physically, such a problem arises, for instance, in the motion of a charged particle in a 
rotating electrostatic potential, as discussed in |12| . 


Brouwer’s particle vs. a particle in a rotating potential. If the graph of the 
saddle surface in Figure[T]is given by a function z = f{x, y, t), then the potential energy 
of the unit point mass on the surface is U(x, y, t) = gf{x, y, t). If we now take the same 
function U as the potential energy of a particle that lives in the plane, we get a problem 
related to Brouwer’s, but not an equivalent one. To note just one aspect of the difference 
between the two problems, note that, unlike a particle in the plane, a particle on the 
surface feels centrifugal velocity-dependent forces due to the constraint to the surface. 
For the motions near the equilibrium these forces are quadratically small (as we will 
explain shortly and as was shown by Brouwer), and they disappear in the linearization. 
And if U itself happens to be quadratic in x and y (as it is in our case, as discussed in 
the next section on the derivation of ([TJ), then the equations for the two problems are 
the same. 
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A particle on the surface vs. a particle in a potential. We would like to add a 
general remark on the difference between two similar problems: a particle on a surface 
z = f{x,y,t) in a constant gravitational field on the one hand, and a particle in the 
plane with a potential U{x,y) = gf{x,y,t) on the other. The two systems have the 
same potential energy - this is the similarity. But here is the difference: the kinetic 
energy for the particle in the planar potential is simply Apiane = while 

2ifsurface = (1 + fl)x^ + + (1 + /^)y^; 

here and fy denote partial derivatives of the height function /, is a much more 
complicated expression. However, since at the equilibrium fx = fy = 0, these derivatives 
are small near an equilibrium, and therefore A^surface ~ Apiane there. Thus the near- 
equlibrium motions of the two systems are quite similar. 

Precessional motion in the rotating saddle trap. It has been known since 
Brouwer that the motion of the particle on the saddle is stabilized for all sufficiently 
high oj [ISIIIT]. An illustration of a similar (but not equivalent) counterintuitive effect 
consists of a ball placed on a saddle surface rotating around the vertical axis and being 
in stable equilibrium at the saddle point HU. We note however that the rolling ball 
on a surface is a non-holonomic system, entirely different from a particle sliding on a 
surface [251 [3U. 

The particle trapped in the rotating saddle exhibits a prograde precession in the 
laboratory frame as illustrated in Figure El This means that the particle moves along 
an elongated trajectory that in itself slowly rotates in the laboratory frame with the 
angular velocity LOp in the same sense as the saddle. Up to now this precession has 
been explained by analyzing explicit solutions of the linearized equations (dniiiuiii]), 
leaving the underlying cause of this precession somewhat mysterious. 



Figure 2: Prograde precession of a particle on a rotating saddle in the non-rotating 
frame for large angular velocity cu calculated according to equations dU). In 
this illustration cj = 20/9 ~ 2.2 (not even that large!). 
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Deriving the equations of motion. 

The motion of a point unit mass x = (x,y) in any time-dependent potential U{x,y,t) 
is given by 

i=-Vt/(x,t) (4) 

(the gradient here is taken with respect to the (x, y)-coordinates). The graph of our 
U (x, t) is obtained by rotating the graph of Uo{x, y) = \{x‘^ — y^), and we must (i) find 
the expression for this rotated potential U{x,t), and (ii) compute VC/ - hopefully in an 
elegant way, without brute force. 

The answer to (i) is simply 


U{x,t) = Uo(R-^x), 


(5) 


where R = R{ujt) denotes the rotation through the angle ojt around the origin (coun¬ 
terclockwise if t > 0) and is given by the matrix 


R = R{iot) 


cos cut — sin ujt 
sin ujt cos iot 


( 6 ) 


Figure [3] explains why: Imagine looking down upon the graph of U{x,t), Figure [3](A), 



focusing on a point x; let us now rotate the graph together with the point x clockwise 
by angle wt, thereby turning the graph into the graph of Uq] the result is shown in 
(B). Since the graph and the point x rotated together, the height above x remained 
unchanged - exactly what ([5]) says! 

(ii) to hnd VC/ we could write ([5]) in terms of x and y and then differentiate, but 
here is a more elegant mess-avoiding way: let us write C/o(x) = ^(x^ ~ y^) as the dot 
product: 

t^o(x) =-(yx,x), where y = diag(l,-l) (7) 

(note that p is the mirror reflection in the y-axis), and use this in jS]): 

Uixi,t) = ^{pR-^x, R-^x) = ^{RpR-^x,x) = ^(S'x,x); (8) 

in dH]) we used the fact that R is an orthogonal matrix. Multiplying out RpR~^ = S 
gives the matrix S defined by equation ([3]). 
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Figure 4: (A): S' is a reflection followed by a rotation. (B). Another interpretation of 
the rotating saddle: rotating the potential with angular velocity u amounts 
to rotating each arrow of the vector field {x,—y) with angular velocity 2uj. 

Now for any symmetric matrix S, one has V(Sx, x) = 2Sx, as one can see with 



V[/(x, t) = Sx; 


(9) 


Equation of motion Q turns into Q, as claimed. 

A “spinning arrows” interpretation of rotating saddle. Note that S{(jjt) is a 
composition of the reflection with respect to the x-axis and the counterclockwise ro¬ 
tation through angle 2ujt, Figure lUj^A). Therefore, for a fixed x and increasing t, the 
vector S{u}t)x rotates counterclockwise with angular velocity 2a;. This leads to the 
following nice geometrical interpretation of the governing equations o-©. The force 
held —S(a;t)x in our equations ([2]) can be thought of in the following way, Figure IHB). 
Starting with the stationary saddle vector held {x,—y), we rotate each arrow of this 
held with angular velocity 2 a; counterclockwise; the result is our time-dependent vector 
held — VC/(x, t) = —S{ujt)x. 

Applications, connections to other systems. 

Before getting to the point of this note, we mention that equations arise in numerous 
applications across many seemingly unrelated branches of classical and modern physics 
[in da da IS]; here is a partial list. They describe stability of a mass mounted on 
a non-circular weightless rotating shaft subject to a constant axial compression force 
da da, in plasma physics they appear in the modeling of a stellatron - a high-current 
betatron with stellarator helds used for accelerating electron beams in helical quadrupole 
magnetic helds |10[ fT^ I18| . In quantum optics, equations ([I]) originate in the theory 
of rotating radio-frequency ion traps m- In celestial mechanics the rotating saddle 
equations describe linear stability of the triangular Lagrange libration points L 4 and L 5 
in the restricted circular three-body problem dsiiiQiisi]. In atomic physics the stable 

^Indeed, the definition of the gradient of a function / states: (V/(x), v) + sv)s=o for 

all vectors v. In our case, /(x) = (5x,x) and -^fix + sv)s=o = (2S'x, v); comparing this with the 
left-hand side of the definition proves \7{Sx,x) = 2Sx. 
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Lagrange points were produced in the Schrodinger-Lorentz atomic electron problem by 
applying a circularly polarized microwave field rotating in synchrony with an electron 
wave packet in a Rydberg atom |21) . This has led to a first observation of a non¬ 
dispersing Bohr wave packet localized near the Lagrange point while circling the atomic 
nucleus indefinitely |22| . Recently, the rotating saddle equations © reappeared in 
the study of confinement of massless Dirac particles, e.g. electrons in graphene [23]. 
Interestingly, stability of a rotating flow of a perfectly conducting ideal fluid in an 
azimuthal magnetic field possesses a mechanical analogy with the stability of Lagrange 
triangular equilibria and, consequently, with the gyroscopic stabilization on a rotating 
saddle |24| . 


The result: a “hodograph” transformation 

The main point of this note is to show that the rapid rotation of the saddle gives rise 
to an unexpected Coriolis-like or magnetic-like force in the laboratory frame] it is this 
force that is responsible for prograde precession. To our knowledge this is the first 
example where the Coriolis-like force arises in the inertial frame. 

To state the result we assign, to any motion x = x(t) satisfying ([2]), its “guiding 
center”, or “hodograph’H 

u = x— —S{ut){'x. — eJic), (10) 

where e = and where J is the counterclockwise rotation by 7r/2: 

■'=(? "J)' 


We discovered that for large cj, i.e. small e, the guiding center has a very simple 
dynamics: 


3 2 

ii - ^ Jii + = 0(6"^); 

4 4 


( 12 ) 


that is, ignoring the e^-terms, u behaves as a particle with a Hookean restoring force 
— and subject to the Coriolis- or magnetic- like force ^Ju. Figure [5] shows a 
typical trajectory of the truncated equation 


ii-Jii H-u = 0. 

4 4 


(13) 


This is in fact the motion of a Foucault pendulum, and just as in the Foucault pendulum 
1261 the Coriolis-like term is responsible for prograde precession of u and thus of its 
“follower” X. 

®The conventional hodograph transformation involves the derivatives of the unknown function. 
Hamilton referred to the path of the tip of the velocity vector of an orbiting planet as the hodograph 
(the word is a combination of the Greek words for path and for trace, or describe.). Following Hamilton, 
in meteorology the hodograph is the trajectory of the tip of the velocity vector; in partial differential 
equations, Legendre’s hodograph transformation involves partial derivatives. We refer to (nni as a 
hodograph transformation since u is a certain mixture of position and velocity. 
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Angular velocity of precession in (IlSp compared to earlier results. Angular 
velocity of precession of solutions of the Foucault-type equation (fT^ turns out to be 
ojp = e^/ 8 . Indeed, writing the equation in the frame rotating with angular velocity 
Wp gives rise to a Coriolis term and a centrifugal term, and the system in that frame 
becomes 

gS ^2 

z + 2ojpJz——Jz + —z + uipZ = Q. (14) 

With the above value of Wp, the second and the third terms cancel, and the resulting 
system 

z + (e^/4 + e®/64)z = 0 (15) 

exhibits no precession (all solutions are simply ellipses). We conclude that e ^/8 is 
indeed the angular velocity of precession of u. This simple expression fits with the 
earlier results and also gives the precession speed for the near-equilibrium motions of 
Brouwer’s particle, as we show next. 



Figure 5: The motion of the “guiding center”, governed by (fT3l) . is the same as that of 
a Foucault pendulum. Here e = 0.45. 


The equations of a rotating radio-frequency ion trap obtained in m reduce to our 
equations (JT]) with e = (jJ~^ = \/ 2 g after the re-scaling of time: t = r-\/ 2 g, where t is the 
dimensionless time in ([T]) and r and q represent here, respectively, the dimensionless time 
and a dimensionless parameter of the trap in |12| . Calculating the angle of precession 
LJpt, we find 

^3 42 2 

Wpt =( 16 ) 

which yields the precession rate obtained by Hasegawa and Bollinger |12) : 




(17) 


Similarly, the dimensional precession frequency of the particle in Figure [T] is: 


Hp — 




(18) 
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3 

A Coriolis—like force in the inertial frame. If were a true Coriolis force, 

it would have been due to the rotation of the reference frame with angular velocity 

3 3 

Y ~ but our frame is inertial. Alternatively, one can think of ^Ju as the Lorentz 
force exerted on a charge (of unit mass and of unit charge) in constant magnetic field 

3 

B = ^ perpendicular to the plane. Rapid rotation of the saddle gives rise to a virtual 
pseudo-magnetic force (cf. |27[ [28112^ 1: as one implication, the asteroids in the vicinity 
of Lagrange triangular equilibria behave like charged particles in a weak magnetic field, 
from the inertial observer’s point of view. 



Figure 6: A trajectory of the guiding center u governed by (|13p (in red) tracking the 
corresponding trajectory x governed by ([2|). The view is in the inertial frame, 
e = 0.45. 


A numerical illustration. Figure [6] shows a solution x alongside its “guiding center” 
u. We note, as a side remark, that near the origin the solution follows the trajectory 
of the guiding center rather closely, reflecting the fact that oscillatory micromotion is 
small near the origin, as is clear from the governing equations ©• 

We discovered the hodograph transformation (IlOp via a somewhat lengthy normal 
form argument mu which, due to its length, will be presented elsewhere. Nevertheless, 
once the transformation (jlOp has been produced, the statement (I12p can be verified 
directly by substituting (1101) into ©; we omit the routine but slightly lengthy details, 
but give a geometrical view of this transformation. 

A geometrical view of the hodograph transformation. Our result says, in effect, 
that the “jiggle” term 

£2 

—5'(a;t)(x — eJx) (19) 

in m, when subtracted from x, leaves a smooth motion0 Figure [7] gives a geometrical 
view of this term. It is still an open problem to find a simple heuristic explanation 

^as far as the powers up to e® are concerned. 
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of choice of the term —eJic in (1101) . Finding a heuristic explanation of the “magnetic” 

3 

term ^Ju remains an open problem as well. 



Figure 7: A geometrical representation of the “jiggle” term S{ujt){x — eJx) (without 
the e^/4 factor) in the hodograph transformation. One should “read” the figure in the 
following order: x x —eJx x — eJx i— S{x — eJx). A geometrical interpretation 
of S is shown in Figure |4l(A). 


Validity of the truncation. Our result says that the fictitious particle u, which 
shadows the solution x, is subject to two forces — (e^/4)u — (e^/4)u, plus a smaller 
0 (e^)-force. 

What is the cumulative effect of this force? One can show, using standard results of 
perturbation theory, that neglecting the 0(e^)-term causes the deviation less than cie^ 
over the time |t| < C 2 £~‘^ for some constants ci, C 2 , for all e sufficiently small. As it often 
happens with rigorous results, this one is overly pessimistic: computer simulations show 
that “sufficiently small” is actually not that small: for example, e = 0.45 in Figure [6l In 
fact, the reason for such an unexpectedly good agreement is the fact that the error on the 
right-hand side is actually 0(e®) - two orders better than claimed, as follows from an 
explicit computation by Michael Berry, m- We do not focus on the analysis of higher 
powers of e because it would only add higher-order corrections to the coefficients on the 
left-hand side of (jl2l) . without affecting our main point (namely, that an unexpected 
Coriolis-like force appears) H 


A proposed experiment. 


Figure [8] illustrates a possible mechanical implementation of the rotating saddle trap (cf. 
[mill])- As we had mentioned in the introduction, a ball rolling on the rotating saddle 
surface m is not the right physical realization of the rotating saddle trap because 

®According to M. Berry’s computation m based on explicit solution of ©, replacing the coefficient 
e^/4 in (ITHll by 


£ 


2 


4 


^ Aa ^ ^ V (1 - VT^){1 - yTTA) 

16 128 J £2 


( 20 ) 


and replacing the coefficient £®/4 by 


£ 

T 


^+ 16 " + 


21 


£“ + •■•=-- 


2 VI - A + VI + A 


( 21 ) 


yields the equation for the exact guiding center (the expression for which is then given by (IIOII which 
includes higher order terms). 
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Figure 8: A possible mechanical realization of the rotating saddle trap. Here x, y are 
the angular variables, and the graph of the potential energy in terms of the 
angular variables x, y is shown. 


(i) the friction is very hard to eliminate, and, perhaps more importantly, because the 
rolling ball does not behave as a sliding particle. In fact, the rolling ball can be stable 
even on top of a sphere rotating around its vertical diameter m'- 

A light rod with a massive ball in Figure [ 8 ] is essentially an inverted spherical 
pendulum; the sharpened end of the rod, resting on the center of the platform, acts as 
a ball joint with the horizontal plane. The springs are attached to the rodjl The height 
of the ball is adjustable, like in a metronome. If the ball is placed sufficiently low then 
the two springs will stabilize the pendulum in the x-direction. Figure [ 1 B )0 and the 
potential acquires a saddle shape (since the y-direction is unstable): 

Uo{x,y) = ^ (ax^ - by^), ( 22 ) 

with a, b > 0, ignoring higher powers of x, y. Here b = g/L, where L is the distance 
from the ball to the sharpened end of the rod. In the next paragraph we suggest a 
simple way to adjust L to make the two curvatures equal: a = b, so that the linearized 
equations become 

x + 65(nt)x = 0, b = Y (23) 

IJ 

a rescaled version of (HD- By rescaling the physical time t to the dimensionless time 
r = ^/at, and introducing 

uj = Q./y/a, (24) 

we transform the above equation into the dimensionless form ([ 1 ]) (after renaming r back 
to t), which is stable if cj > 1 [Tj [iSl |T7]. Now a = b = gjL, where L is the distance 
from the mass to the ball joint, and (|24l) gives us the length 

L = (25) 

®Theoretically, we want to avoid transferring the rotation of the platform to the rotation of the rod 
around its long axis (thereby affecting its dynamics), and thus the attachment should, theoretically, be 
via some frictionless sleeve. In practice, however, this frictionless sleeve will hardly make a difference: 
the moment of inertia of the rod+ball around the long axis is negligible compared to the moment of 
inertia relative to a diameter of the platform, and thus the dynamics of the rod will be little affected 
by its axial rotation. 

^here x and y are angular variables. 
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The 78 rpm of a vinyl record player corresponds to 0 ~ 8.2sec“^, and the value oj ~ 2.2 
referred to in Figure [2] corresponds to the height L ~ 71 cm. How short can we make the 
pendulum without losing stability? The cutoff length is L ~ 14 cm, as follows from (|25p 
and the fact that ([T]) is stable if and only if w > 1. Incidentally, large L corresponds to 
large w, according to (f25ll (for fixed rpm). This makes intuitive sense, since a “natural” 
unit of time in our system is the period 2'K-\/L/g of the oscillations along the stable axis 
of the stationary potential; for large L this period is long, and during it the rotating 
potential will spin many times, corresponding to large oj. 



Figure 9: How to find the length for which the curvatures of the saddle are equal and 
opposite: Uxx{0,0) = —Uyy{0,0): hnd by trial and error the height He at 
which the x-direction changes stability, and half He- 


How to (easily) realize the saddle with equal principal curvatures. It may 

seem (as it did to us initially) that one needs to measure the stiffnesses of the springs, 
the various lengths in Figure IH the mass of the ball, and then use these to compute the 
value of L. Instead, here is a way to avoid all this work. Referring to Figure [9l adjust 
the position of the massive ball along the rod to such critical height He as to make the 
ball neutrally stable in the x-direction: if the ball is too high, it will be unstable in 
the x-direction; if the ball is too low, it will be stable in the x-direction; a bisection 
method will quickly lead to a good approximation to He- Remarkably, the desired 
“equal curvatures” height of the ball is simply 

L = ^. (26) 

An explanation. The angle y satisfies the inverted pendulum equation y — ^ sin y = 
0, which for small angles is well approximated by 

y - fy = 0. (27) 

Similarly, the linearized equation for the angle x is 

X + + ^^x = 0; (28) 

L is not yet chosen and k depends on the parameters of the system, but not on 
Our goal is to hnd L such that the coefficients in the above equations are equal and 

®Specifically, k depends on the springs’ stiffness, on the angle they form with the diameter in Figure|8] 
and on the ball’s mass (to which it is inversely proportional). But an easier way to find k is to note 
that k = gHc, as we show below. 
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( 29 ) 


opposite, which amounts to asking for L satisfying f: = or 

2Lg = k. 

We now relate the unknown k to He- For the pendulum of the length He the angle x 
satisfies x + {--^ + ^)x = 0, with the coefficient + ^ =0 since the equilibrium 
is neutral; this gives 

k = gHc- (30) 

Substituting He = 2L into the last equation gives (I29p - precisely the condition for the 
equality of the coefficients in U. 

Conclusion. We showed that the rapid rotation of the saddle potential creates a weak 
Lorentz-like, or a Coriolis-like force, in addition to an effective stabilizing potential - 
all in the inertial frame. We also proposed a simple experiment to demonstrate the 
phenomenon. 
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